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TECHNICAL MEMORANDUM X- 64746 


SKYLAB ATTITUDE CONTROL AND ANGULAR MOMENTUM 
DESATURATION WITH ONE DOUBLE-G IMBALED 
CONTROL MOMENT GYRO 

INTRODUCTION 


The Skylab is a scientific space station in circular earth orbit with 
more or less cylindrical moment-of-inertia distribution. The minimum 
moment-of-inertia axis is nominally in the orbital plane and perpendicular 
to the sunline. The attitude is normally maintained by three double-gimbaled 
control moment gyros (CMG) [ 1] - Several possibilities exist on how to 
control the attitude when two CMG's have failed and normal attitude control 
and angular momentum desaturation [2, 3] have to be abandoned. The 
possible attitude control concepts are the following: 

Concept A. The nominal attitude reference is retained; the remaining 
CMG is also deactivated and the thruster attitude control system (TACS) 
controls all axes to the specified attitude and rate deadbands. 

Concept B. The nominal attitude reference is retained; the remaining 
CMG controls sun-pointing (two axes) accurately and the TACS controls the 
third axis within a specified attitude and rate deadband. 

Concept C. The attitude reference follows an unstable limit cycle of 
about 0. 3 rad amplitude which the vehicle exhibits if it is allowed to yield 
under the influence of gravity-gradient torques; the remaining CMG is also 
deactivated and the TACS controls all axes to the specified attitude and rate 
deadbands. 

Concept D. The attitude reference is commanded to oscillate in the 
orbital plane with 0. 2 rad amplitude and twice orbital frequency; the remaining 
CMG controls the commanded attitude about orbital north accurately as well as 
the attitude about the vehicle minimum principal moment-of-inertia axis and 
the TACS controls about the projection of the sunline into the orbital plane, 
but to a rate deadband only. 

The advantages and disadvantages of the concepts are discussed in the 
following paragraphs. The numerical values are based on the vehicle data of 
Tables 1, 2, and 3. 



TABLE 1. GENERAL DATA 


CMG Momentum Magnitude 

H = 

3115. 

3 [ N* m* s] 

TACS Minimum Impulse Bit 

MIB 

= 18 

[N*s] 

Momentum Change in x for one MIB 

AH 

X 

= 60 

[N* m* s] 

Momentum Change in y or z for one MIB 

AH = 

= 210 

[N* m* s] 

Total TACS Fuel Available (in Impulse 
Units) 

*tacs 

= 160 000 [N* s] 


TABLE 2. COMMAND/SERVICE MODULE DOCKED TO SKYLAB 


Minimum Moment-of-Inertia 

I = 0.992* 10 6 [ kg* m 2 ] 

X 

Intermediate Moment-of-Inertia 

I = 6. 168* 10 6 [ kg* m 2 ] 
z 

Maximum Moment-of-Inertia 

I = 6.245* 10® [kg* m 2 ] 

Amplitude of Cyclic Momentum 

H « 4500 [N*m* s] 

(Concepts A and B) 

eye 


TABLE 3. COMMAND/SERVICE MODULE NOT DOCKED TO SKYLAB 


Minimum Moment-of-Inertia 

I = 0.893* 10® [kg* m 2 ] 
x 

Intermediate Moment-of-Inertia 

I = 3.813* 10 6 [ kg* m 2 ] 

y 

Maximum Moment-of-Inertia 

I = 3.882-10 6 [kg*m 2 ] 
z 

Amplitude of Cyclic Momentum 

H * 2800 [N* m* s] 

( Concepts A and B) 

eye 













Concepts A and B require no modification of the attitude error genera- 
tion. However, the cyclic angular momentum caused by the gravity-gradient 
torque is so large that the available TACS fuel is exhausted within 3 days 
(assuming 50 000 [N- s] available at the beginning) for Concept A and within 
about 5 days for Concept B; i. e. , the mission would have to be terminated 
very soon. Concept A would not allow many solar experiments, whereas 
Concept B would allow most. Concept C has been called quasi-inertial attitude 
control [4,5] and for an ideal case the fuel consumption has been reduced 
drastically. The fuel consumption is unfortunately very sensitive to the know- 
ledge of the principal moment-of-inertia axes misalignment and to the know- 
ledge of the attitude of the minimum principal moment-of-inertia axis with 
respect to the orbital plane (requiring a functioning star tracker) and increases 
with the addition of other disturbance torques besides the gravity-gradient 
torque. Only the biomedical experiments can be performed for Concepts C 
and D. The fuel consumption for Concept D, however, can be made completely 
insensitive to information on principal moment-of-inertia axes misalignment, 
and disturbance torques and does not require star tracker information since 
no attitude information is used for the TACS command, but only rate informa- ; 
tion. The worst case fuel consumption can be held below 50 [N- s] per orbit 
for all solar elevation angles with respect to the orbital plane and all projected 
disturbance torques. The control system parameters, however, are assumed 
to be reasonably close to their optimum. This compares very favorably with 
the 100 [N* s] per orbit for Concept C, and the 1200 [N- s] per orbit for Con- 
cept Av It should be pointed out again that the low fuel consumption for 
Concept C only applies to the ideal case where the principal moment-of- 
inertia axes misalignments are known exactly and where there are no disturb- 
ance torques other than the gravity gradient torque. A 50 [N- s] per orbit 
figure is equivalent to 66 days operation (assuming 50 000 [N- s] TACS fuel 
available at the beginning) versus less than 3 days operation for 1200 [N- s] 
per orbit. The attitude reference oscillation is required for Concept D so 
that part of the cyclic momentum is stored in the vehicle, reducing the' 
remainder to the point where one CMG can accommodate it and have some 
margin for control. The details of Concept D are set forth in the body of this 
report. „ 


GRAVITY-GRADIENT TORQUES AND ACCUMULATED MOMENTUM 


When the gravity- gradient torque is expressed about vehicle principal 
axis, one has [3] 


3 



( 1 ) 


T 

-gg 


3fi 2 


(I - I ) r r 

z y y z 


(I - I )r r 

X Z Z X 


(I - I ) r r 
y x x y 



where it is assumed that the space vehicle is in a circular orbit with angular 

velocity ; I , I , I are principal moments of inertia; and r , r , r 
x y z x y z 

are the components of the local vertical in the principal coordinate system. 

To bring out the essentials of single CMG control for Skylab, it can be 

assumed that there is no misalignment between geometric and principal axes 

and that I = I =1 or AI=I -I = -(I -I ) which also allows one to 
y z y x x z 

disregard the solar elevation angle (see Appendix A for definition) since 

now any two mutually perpendicular axes, which are also perpendicular to the 
principal x-axis, are principal axes. With this simplification, one has 



Now, one assumes that the vehicle is rigidly controlled about orbital 
y to the angle 0^ and about vehicle x to zero, while it is allowed to have a 

small error <p in orbital z. Then, one has 
zeo 




where n = n - 0 and n is the angle from orbital midnight. Small angle 
tt t yc l 

approximations have been assumed for 0 with c 0 =1 and 

zeo zeo . 

s0 =0 . Insertion of the r components into equation (2) yields 

zeo zeo 


T 

-gg 


= -3fl 2 AI 


V s \t 


Transformation into orbital coordinates and substitution of (n - 0 ) for 

. , . t yc 

7] fct yields 



C 0 

yc 

0 

S0 

yc 


l. 

-0 

zeo 

o' 


0 

T = -352 2 AI 

-ggo 

0 

i 

0 


0 

zeo 

1 

0' 


S2(v t - (t yc )/2 


— S0 
yc 

0 

C0 

yc 


0 

0 

1 


Keo s2(ri t ~ V* 




2,> zeo s( ’ , t - V )C ”t 

82 W 

2 ' # zeo S< ’ , t - V S, t 


( 4 ) 


The orbital y-torque will be discussed first, since it is not a function 

of the small angle 0 . The angle 0 allows the vehicle to yield some 

zeo yc 

under the influence of the gravity-gradient torques and, therefore, absorb 
angular momentum which otherwise would saturate the one operative CMG. 
The following form is chosen for the orbital y-command (A is a constant 
amplitude) 
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V = As2( ’ l t - V 


yd' 


which yields for the orbital y-torque [refer to equation (4)] 


T = --|s2 2 AIs2[7j - As2(r] - 0 )] 

yo 2 t t yd 


yd 


The angle 0 is a desaturation angle. This is shown in Appendix B 

with the results that the momentum accumulated over half an orbit is (H is 
the initial condition) y ° 1C 


H -H . 
yo yoic 


7r/4 

-6ftAI f c2r] c(K s2?7 f )s(K c2r] )dr? 

^ C C C S i l 


(5) 


with = 2Ac2 <p and K g = 2As20^. The actual integration in equation (5) 

is only from 0 to 7r/4 due to utilization of symmetry properties (see Appendix 

B). H (rj X for three different values of 0. , is shown in Figure 1. When 
yo t yd 

0^ is small, equation (5) simplifies to 


H - H . = G ,0 , 

yo yoic yd yd 


( 6 ) 


7r/4 

with G = - 24S2AI A J c 2 (27J t )c(2As2Tj^)d77 . Equation (6) shows that the 
y 0 

accumulated momentum is proportional to 0 , and the desaturation angle 

yd 

effectiveness G is roughly proportional to the amplitude A . 
y u 

The orbital z-torque of equation (4) is 


T zo = -{3S2 2 Als[rj t _As2(77 t -0 )]sT? t }0 


zeo 


6 



-0.02 [rad] 
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Figure 1. CMG orbital y- momentum component (ideal limit cycle). 



This torque has in general the opposite sign of 0 


zeo 


Only very temporarily 


(for 0.^0 and n small) can it have the same sign as 0 .It is, 
v r yd t zeo 

therefore, a restoring torque which drives the minimum principal moment- 
of-inertia axis toward the orbital plane. However, since the quantity in 
braces is not a constant, instability occurs unless damping is provided. Non- 
linear damping is provided by the fact that thrusters are fired to reduce 

0 whenever it exceeds a specific deadband, 
zeo 


excursions of 0 


zeo 


The result is that the angular 
are kept within acceptable limits. 


The orbital x torque of equation (4) is 

T xo - . 

Since 0 is not directly controlled, no control is possible over the accu- 
r zeo 

mulated x momentum. It has to be absorbed by the C MG and thruster firings 
will keep the momentum within acceptable limits. 

ORBITAL y-MOMENTUM DESATURATION 


The commanded orbital oscillation 0 must be exactly in phase with 

yc 

the direction of the gravity gradient, or bias momentum accumulates as 
shown in the section on Gravity-Gradient Torques and Accumulated Momentum. 
On the other hand, this fact can be used to desaturate orbital y- momentum, as 
well as seek the direction of the principal x-axis. The commanded oscillation 
is 


V = As2 <v V 


(?) 


where is a desaturation angle to be updated every half orbit. Any change 


of 0 will also lead to a change in 0 
yd yc 


This causes a maneuver requiring 
momentum from the CMG. The start of the maneuver, i. e. , the checkpoint. 
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should be placed at a point when ample momentum is available, i.e. , away 
from the momentum peaks. The times when rj^ = 7r/4 and 7^ = 57r/4 have 

been selected as the two checkpoints per orbit. A typical momentum curve 
(docked configuration) for no inertia misalignment and no other than gravity- 
gradient torques is shown in Figure 1. 


In general, a momentum ramp AH^ will be superimposed on the 
cyclic momentum. The basic desire for the desaturation angle change A<£ 
is to eliminate the ramp and simultaneously average the momentum peaks: 


yd 


A<p: , = (A <p ) + (A0 ) 

yd yd. p yd r 

. = k (ah + Aii ) 

.vd p r 


with K = -l/G , . While the steady-state effect of <p ' is given by equa- 
yd yd yd 

tion (-6) , a Af , also has other effects. One is the different momentum 

yd 

stored in the vehicle due to the difference in <p : 

. ■ yc 


A? • = 0 (+) - 0 (-) 

yc yc yc 


= 2fiAc2[7? t - (0 +A0 )] - 2flAc2[7? t - 


( 8 ) 


For sufficiently small A0^ d (c2A^ d ~ 1 , s2A0^ d ~ 2A0 d )', one has 


A <P ~ - 4fiAs2(T? - 0 ,)A0 , 
yc t yd 7 yd 



or 


(AH). , = K a . 

A y d 


with = -4J2I^ c . The existing ramp can be established by 


AH = H - H' / 
r yo yo(n-l) 


( 9 ) 


with H' , ,» = H , . + K. , .A <P , where the (n-l) indicates 

yo(n-l) yo(n-l) A(n-l) yd(n-l) 

the values at the last checkpoint. To average the peaks, one uses 


AH = (H + H . )/2 + (0H /0H )AH 
p max min p r r 


( 10 ) 


The effect of A H on the peals can be broken down into three dif- 
ferent parts: 


(3H /3H ) = K + K K + K 
p r r A yd <f> 


Looking back, the peaks should be adjusted to 


H' = H +[(t-t )/(0.5T )] AH 

max max max orb r 


H' . = H . + ((t-t . )/(0.5T ) ] A H 

min min min orb r 


Both effects can be combined into K AH with 

r r 


K = (2t-t -t . )/T , 
r max mm orb 


( 11 ) 
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Looking ahead, the elimination of the ramp will generate a momentum change 
of and will affect the peaks some more by the phase shift and by 

the change of form of the momentum curve. These effects are lumped into 

KM . Since the value of K depends on the angle 0 , , on the vent 
<p v 0 yd 

torques, and on the value of the ramp, an acceptable K has to be established 

by simulation. 

The final result is 


A0 , = K { (H +H . )/2+ ( 1+K +K K +Kj(H -H' , .)}' 

yd yd max min r A yd 0 y o yo(n-l) 


( 12 ) 


To avoid excessive maneuver momentum requirements or maneuver 
time, A'Pyd is limited to A0^^ • The larger the nongravitational disturb- 
ance torque, the larger the desaturation angle and the larger the peak-to-peak 
distance of the cyclic component of the orbital y- momentum. But the peak-to- 

peak distance (H -H . ) has to be kept below twice the orbital y- 
max min 

momentum limit to avoid TACS firings which do not oppose the disturbance 
torques. This can be achieved by a limit placed on 0 ^ . The value of the 

limit cannot be a constant, since it depends on the inertia misalignment and, 

also, on the disturbance torque level. Therefore, the value of the limit has 

to be changed, when the changing conditions demand it. When the peak-to- 

peak distance is below its allowed maximum H and the desaturation angle 

P p 

0 ^ was not on the limit for the past half orbit, no adjustment of the limit 
0 is necessary. If, however, 0^ was on its limit, the limiting value 

is too low and has to be increased. If the peak-to-peak distance is more than 

H , the limiting value has to be decreased, whether 0 , was limited or 
PP yd 

not. A minimum value is also imposed on the limit itself, since its lowest 

reasonable value can easily be established, but transients could let the limit 

drop lower. 


<P 


ydl 


After 0 , is updated by A 0 , , it is limited to the latest value of 
yd yd 

Details of the logic and equations are given in Appendix C. 
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ATTITUDE ERROR GENERATION AND CMG CONTROL LAW 


A double-gimbaled CMG with fixed angular momentum magnitude H 
has only two degrees of freedom and, therefore, can control only two axes. 
It was decided to select the vehicle geometric reference axis x^ and the 

orbital north axis y^ (see Appendix A for coordinate system definitions). 

Several other x-axis choices could have been made (x , x ) but these 

pr or 

would have resulted in more complex equations and would have required 
additional information and angles (moment-of- inertia misalignment, princi- 
pal Zp-axis elevation). The orbital north axis is determined under the 

assumption that the vehicle x^-axis is in the orbital plane. The resulting 

error is small and can be tolerated. Otherwise, accurate angle information 
about the sunline would be necessary, requiring the star tracker to be oper- 
ative. 

The commanded oscillation about the orbital north is achieved by 
generation of the attitude error and rate error as 


<t> 

xe 


%al" 


‘o' 

0 

ye 

= 2S ^<V4> 

v 2 

- 

e 

y 

0 

ze 


%a3 


e 

z 


> 

xe 


\ — 

i 


’ 0 

ye 

= 


- 

e 

y 

'<p 

ze 


1 

N 

l 


e 

z 


(13) 


( 14 ) 



with 



0 crj 
yc 


x 


£ = _0 S77 

z yc x 


€ = 0 cr) 

. y yc 'x 


e = -0 sn 
z yc x 


and 


0 = As2(n - 0 ) 

yc 't yd 


(15) 


. 0 = 2S2Ac2(tj - 0 ) 

yc t yd' 


(16) 


where (i = 1, 2, 3, 4) are quaternions relating CS to CS , 

A is the amplitude of the commanded oscillation 0 c , SI is the average orbital 
rate, and 0 ^ is an orbital y- momentum desaturation angle (see chapter on 
desaturation) . 


The attitude and rate errors have to be modified as follows: 
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The orbital z-component in 0 and 0 have been suppressed. The 

attitude error can now be limited to a relatively small value, without the 
possibility that a possibly large orbital z-error becomes bothersome (would 
limit the commandable velocity when rj + 0) . 

X 


The normalized CMG torque command can now be generated as usual 

15]-; ' 
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( 23 ) 



(The normally applied bending filter is neglected here. ) 
A CMG steering law of the form 


o> e e 

il il cl 



(known as the cross-product steering law) will provide a torque component 
in the desired direction, but its magnitude will always be smaller than 
desired. This gain reduction can be ignored because the gain drops only to 
70 percent (TACS firings prevent lower values) for the worst case orbital 
desaturation transient and is normally higher than 95 percent. While the CMG 
will provide x- and y-control, the necessary change of its z-momentum 
component will act as a disturbance on the z-axis. 

The cyclic orbital y-momentum and the commanded oscillation about 
orbital y (it results in an oscillation of the CMG angular momentum vector 
in vehicle space) can be accommodated by a "window” of about ± tt / 12 width in 
the orbital plane and about ±7r/4 height along orbital north. This window is 
enforced through TACS firings (see TACS control). The center of the window 
has been selected as follows: (for T] x = 0) 

For CMG No. 1 operative 
+ 0.479 
0 

-0.877 

6 l(l) = + ( 7r “ 1 ) /2 = L071 [rad l 



(25) 
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For CMG No. 2 operative 


1 

CD 

ro 


l 

O 

1 

e 22 

= 

0 

e 23 


1 j 

H- 4 . 


6 1 ( 2 ) - ° 

6 3(2) = 1/2 [radl 


( 26 ) 


For CMG No. 3 operative 


e 31 


0.479' 

e 32 

. = 

o' 

e 33 


-0.877 


6 1(3) - ° 

6 o fr>\ = + ( 7r - 0-5) [rad] 


(27) 


This window selection insures that no gimbal stops are encountered for the 
whole range of possible solar elevation angles 77 (-1.265 [rad] to +1.265 

[rad]). 


TACS CONTROL 


The thruster attitude control system (TACS) is used to insure that the 
CMG angular momentum does not exceed the prescribed "window" and it also 
insures that the vehicle angular rate about the axis not controlled by the GMG 
remains within a specified deadband. 
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The TACS consists of six cold gas thrusters, which are selected and 
fired individually. The thrusters are mounted on Skylab as shown in Figure 2. 
The main torque polarities of the thrusters are shown in Table 4. It is also 
shown in the table that there is a built-in cross -coupling between the vehicle 
x- and z-axis, and only firing of the appropriate pairs will eliminate this 
cross-coupling, as much as tolerances and misalignments allow. 

The thruster attitude control system fires a minimum impulse bit 
whenever ||i,| s 1 for any axis [ 6 1 - This property is used in Table 5 by 

allowing only three states for 0 i 0 and ±1. 

The thrusters are selected according to Table 5. 

For proper initial conditions and for no other than gravity-gradient 

torques, the cyclic gravity-gradient torques are absorbed by the CMG and 

the nominal momentum profile of Figure 1 results. The variation of the 

angular momentum of the CMG along the orbital z-axis (see Figure 3) acts as 

a disturbance and the gravity-gradient torque provides a restoring torque on 

the minimum moment-of-inertia axis whenever a z-torque drives it out of the 

orbital plane by 0 . The nominal profiles for <p and <p of Figures 

zeo zeo zeo 

4 and 5 result if the proper initial conditions have been chosen. The nominal 

0 profile is an unstable limit cycle and will diverge until 0 exceeds 
zeo zeo 

its deadband and is reduced by a thruster firing (Figs. 6 and 7). The objec- 
tive of minimum fuel consumption requires that the deadband is wide enough 
to allow the ideal limit cycle, but still narrow enough such that only one firing 
is required whenever TACS is actuated, which should reset the initial condi- 
tions as close as possible to the ideal case with the hope that there will be 
several orbits devoid of thruster firings. In other words, the value of the 
deadband is very critical and will have to be set repeatedly during flight, 
especially since it is also dependent upon the momentum change due to a 
minimum impulse bit from the thrusters, as well as disturbance torques and 
other influences. While a wrong deadband value (either too high or too low) 
is not disastrous, it does result in an increase of the thruster fuel consump- 
tion, and should be corrected as soon as possible. 1 

A constant z-torque results automatically in an average hang-off of 
the x-axis out of the orbital plane. Several times the presently assumed vent 
torques can be handled without an increase in the fuel consumption. The 
deadband, however, has to be adjusted if the torque varies. 


1. The worst-case fuel consumption of 50 [N-s] per orbit (quoted in the 
introduction) does not assume worst case for the deadband, but only for the 
solar elevation angle and the disturbance torques. 
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Figure 2. TACS configuration. 






TABLE 4. MAIN TORQUE POLARITY 


Thruster 

x -Axis 

V 

y -Axis 
J v 

z -Axis 

V 

1 

0 

- 

0 

2 

- 

0 

- 

3 

■ + 

0 

+ ' 

4 

0 

+ 

0 



f\ 


<3 




6 


0 

+ 


TABLE 5. THRUSTER SELECTION 


Condition 

Selected Thruster 


+ 1 

2 

i l> ' 

0 

None 

X 




- 1 

3 


+ 1 

1 


0 

None 

y 




- 1 

4 


+ 1 

5 

4> 

0 

None 

z 




- 1 

6 


In the following discussion, reference to an x-firing indicates the 
firing is caused by the momentum exceeding the limits set for , a y-firing 

is caused by exceeding its limits, and a z-firing is caused by ^ zeo 

exceeding its deadband. 

A constant y-torque results in the appropriate change in the desatura- 
tion angle <p and in an increase in the amplitude of the cyclic momentum, 
yd 
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Figure 4. Angular velocity ^ zeQ (ideal limit cycle). 


0.01 



Figure 5. Angle £ (ideal limit cycle). 



0 000 1 
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Figure 6. Long term instability of 



Since the amplitude is bounded, the desaturation angle <^> has to be limited, 
and any excess momentum has to be desaturated by the thrusters. However, 
the <2>yd-limit can desaturate about twice the anticipated disturbance momentum. 


A constant x-torque will have to be eliminated entirely by thruster 

firings. Since the z-limit cycle is very sensitive, any cross-coupling from 

x-firings would be highly undesirable and therefore pure x-firings will always 

be in pairs. However, when a thruster is fired on a d> command, it is also 

z 

utilized to reduce the x-momentum with respect to the nominal x-momentum. 

Due to the commanded oscillation in the orbital plane, the nominal x- 

momentum translates into a sinusoidal momentum H , in vehicle coordinates: 

xb 



Hs(<j!> - <p 

yc 



(28) 


where <p- ^ is a bias angle that depends upon which of the CMG’s is operative 

(<p ■■ = -0. 5 [ rad] for CMG 1; <p . = 0 for CMG 2; 0 , = +0. 5 [ rad] for CMG 
yb yb yb 

3). 


A y-firing (indicated by i ip * 0) will, through its effect on the CMG 

orbital z^momentum, also always decrease the magnitude of <p .It should, 
• > zeo 

therefore, receive priority over the other firings. If I?? I > tt/4 , the appro- 

X 

priate polarity ip is selected instead of ip . A z-firing is not only called 
z y 

for whenever <p exceeds its bounds, but also when it can be clearly estab- 
zeo 

lished that it has the wrong polarity with respect to the desired limit cycle. 

Figure 3 in conjunction with Figure 1 shows that <f> should be positive for 

zeo 

Ih I less than some lower value H ,, and negative for |H I above some 
’ yo 1 yll & yo 1 

upper value H The values for these limits depend on the disturbance 

torques and have to be established during flight. Once the need for a z-firing 
has been established, it remains to be determined which conbination of ip 

y 

and ip (which produce momentum changes in vehicle coordinates) is most 
z 

beneficial in reducing the magnitude of <j> while minimizing the cumulative 

zeo 

effect of the z-firings on H . The reason for the latter is elimination, as 

yo 
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much as possible, of cross-coupling into . To do this, five firing 

combinations are evaluated in terms of their effect on H and b . The 
' , . . yo zeo 

five combinations are 



0 

+ 1 

- 1 

+ 1 

- 1 


i p 

r zo 

ip 

T zo 

Ip 

T zo 

0 

0 


The following formulas are used (the latter is also in momentum rather than 

(29) 

(30) 

(31) 


H 2 - H 2 - (H + AH ) 2 
xo yo yo 

- / H 2 - H 2 Hi 2 . (32) 

n/ x o yo 

Equation (29) represents the angular minimum impulse bit for a vehicle y- 

or z-firing. R is the thruster lever arm (a negative value) , f is the 
X c 

thruster force. At is the commanded ON- time and t., is a correction 
mib td 

factor for thrust decay. AH' gives the nominal AH for a single engine 

zo 

firing and it is used in equation (33). The difference in the roots of equation 



angular velocity) : 




AH = R f (At + t. ) 
x t mib td' 


AH' = AH[max(cr? , | S 77 |)] 

• X X 


AH = AH(ct? ip - stj ib ) 
-yo x y x z' 


AH = AHc0 (sti ib + cri ib ) 
zo yc x y x z 
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through 


(32) represents the effect of the orbital y-momentum change on i > zeQ 

the CMG momentum change. Any firing combination which does not produce a 

AH with at least half of AH’ with the correct polarity will be disregarded, 
zo 

The rest will be evaluated according to their demerits; 


Yah + a h 

+ ^ 

AH tp - AH' 

LJ yoz yo 

d 

zo zo 


(33) 


- P 

where >AH is the accumulated effect of the z-firings on H since the 
yoz yo 

last checkpoint. The first term results in demerits for increasing the effect 


on H 


yo 


the last term penalizes deviation from the most beneficial single 


engine firing, and (<l) allows relative weighing of the two. The firing 
combination with the smallest demerit D is selected. 


Simpler firing logic has been investigated, but large cross-coupling 
was produced between the orbital y- and z-axis, resulting in strong positive 


feedback at certain 
teed. 


7 ]^ angles, so that settling out of <f> was not guaran- 


A flow chart of the logic with equations is given in Appendix D. To 
reduce coupling into the desaturation equations, the following quantities have 
to be updated by the expected value of AH (= is here a FORTRAN equal): 


H' 1 = H’ , > + AH 

yo(n-l) yo(n-l) yo 


H 

H + AH 


max 

max yo 


H . = 

H . + AH 


min 

min yo 


Yah 

= Ya H + 

AH 

^ yoz 

4 -> yoz 

y° 

Yah 

= Yah + 

AH 

u yoy 

U yoy 

yo 
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DISCUSSION 


Long-term drift of the strapdown computation for the attitude reference 
can be eliminated by slightly modifying the normal Skylab updating method, 
where the output of the sun sensors is compared with the calculated attitude. 
Updating is possible whenever the magnitude of the sun sensor output is below 
0. 1 rad. This occurs three times per orbit (sunrise terminator, noon, and 
sunset terminator) and lasts for about 460 s (for an amplitude A of 0. 2 rad). 

Some TACS ONLY attitude control along the lines of References 4 and 
5 must be incorporated to be prepared in case the last of the three ATM CMG’s 
also fails. 



APPENDIX A 


DEFINITIONS OF SKYLAB COORDINATE SYSTEMS (OS's) 


Only the coordinate systems (OS’s) needed for this memorandum are 
defined. All OS’s are right-handed and orthogonal. All angles are defined 
mathematically positive. • 


QvrmKnl 

*w»jr l 


ni Hnn 


X 

o 


N/A 


Basic orbital CS; z toward 
o 

ascending node; y toward 
orbital north. 


z 

o 


y 


x = [77 ] [77 ]x 
r x y o 


Reference CS; toward sun; 

x in the orbital plane; y in 
r r 

northern orbital hemisphere. 


or 


or 


or 


X = [77 px 
or x r 


Orbital references CS; z along 

or 

projection of z^ into orbital 

plane; y toward orbital north, 
or 


z 

v 


X = [q]X 
v Va r 


Vehicle geometric CS. Also the 
attitude control system CS. 

[ 0 va 3 indicates the attitude 

deviation from x (due to error 
r 

or commands). This CS is also 
assumed to be the principal CS 
for the development of the equa- 
tions in this report. 
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Symbol 

Transformation Matrix 

Definition 

X 

VC 



Commanded vehicle CS. n in 

X 

y vc 

X 

VC 

[r? 1 [ 0 ]'[0 IX 

x zeo yc or 

this transformation is only 
numerically equal to the solar 




elevation angle, but physically it 

Zi 

VC 



is a different angle. : 

x d 



Disturbance CS; z , toward center 
d 

y d 

x d - 

[7?JX 
t or 

of earth; y^ toward orbital north 
( r = 0 indicates orbital mid- 

Z d 



night; rj is ay angle). 
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APPENDIX B 


DESATURATION EFFECTIVENESS 


The gravitational orbital y-torque is [refer to equation (4)] 


T = --S2 2 AIs2[rj t - As2(r? i - <p J] 


yo 2 




where is the desaturation angle. The y- momentum after half an orbit 

is (H . is the initial momentum) 
yoic 


T . /2 

orb 


H = H . + / T dt 
yo yoic J yo 


or 


IT 


H - H = - — fiAI / s2[n - As2(n - <p ,)]drj 
yo yoic 2 0 J 't w t yd' 't 


The integral is broken into two parts 


H - H . = - — £2AI (I t + I 2 ) 

yo yoic 2 


with 


■n / 2 

h = J s2[rj t - As2(7j fc - 0 yd )]dT? t 


n 

I 2 = y s2[TJ t - As2(r? t - 0 yd )]d77 t 
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A substitution of a - n - 


into I 2 yields 


tt/2 

1 2 = ~ J s2[o! - As2(a + <p )]cta 

n y u 


Since a is only a dummy integration variable, any variable can be substituted. 
The author elected to substitute 77 and get for the sum of the two integrals 

7r/2 

(H y " H yo )/( - 3S2AI/2) = J F ^ t )dT] t 

with 

■F.(T? t ). = s2[i? t - As2(r? t - 0 yd )I - s2[r/ t - As2(?7 t + 0 )] 


This expression shows that for 
development of F^) yields 



the integral is zero. 


Further 


F(?7 t ) = s2{r} t -A[s2rj t c20 yd -c2rj t s20 yd ]} - s2{77 t -A[s27J t c20 yd +c2rj t s2^ d ]} 


= s2(i7 t -As2r? t c2^ d )c(2Ac27} t s20^ d ) + c2(rj t -As2T7 t c20 yd ) s(2Ac2r) t s2<£ yd ) 


-S2(rj -As2rj c20 )c(2Ac2t? s2<p ) + c2(r? -As2rj c2^ ,)s(2Ac2rj s20 ,) 
t l ya t ya t t yd t yd 7 

= 2c2(t] -As2ijc2</> ) s(2Ac27? s2<p .) 

t t yd t yd' 

= 2[c2rj c.(2As2rj c20 ,) s( 2Ac2rj,s20 ' 
t t yd t yd 

+ s2tj s(2As2tj c20 )s(2Ac2r] s2<£ )] 
t t yd t yd 


32 



The last tern of F(t] ) is antisymmetric with respect to r\ = 7r/4 and does 

V C 

not result in a contribution when the integration is performed from = 0 to 
t^ = 7t/ 2. The final result is 

rr/4 

c2t) c(K s2t),)s(K c2r) )dn 
L c c s t t 

with 


H - H . = -6QAI 

yo yoic 


K = 2Ac20 . 
c yd 


K = 2As2 <p , 
s yd 

where the integration is only from zero to tt/ 4 because of the symmetry with 
respect to tt/4 and the result is doubled. The effectiveness of 0 had to 

be established by numerical integration, but small angle approximation shows 
that the accumulated momentum is proportional to 0 : 


H - H . = G 0 

yo yoic yd yd 


with 



it/ 4 

-24S2AIA J c 2 (2t? )c( 2 As 2?7 )dr; 

0 1 1 
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APPENDIX C 


LOGIC FOR ORB ITAL y-MOMENTUM DESATURATION 
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APPENDIX D 

EXTREMA OF H AND TACS FIRING LOGIC 
yo 
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